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ABSTRACT 

A forward model is described in which we synthesize spectra from an ab-initio 3D MHD simulation 
of an outer stellar atmosphere, where the coronal heating is based on braiding of magnetic flux due to 
photospheric footpoint motions. We discuss the validity of assumptions such as ionization equilibrium 
and investigate the applicability of diagnostics like the differential emission measure inversion. We 
find that the general appearance of the synthesized corona is similar to the solar corona and that, on 
a statistical basis, integral quantities such as average Doppler shifts or differential emission measures 
are reproduced remarkably well. The persistent redshifts in the transition region, which have puzzled 
theorists since their discovery, are explained by this model as caused by the flows induced by the 
heating through braiding of magnetic flux. While the model corona is only slowly evolving in intensity, 
as is observed, the amount of structure and variability in Doppler shift is very large. This emphasizes 
the need for fast coronal spectroscopy, as the dynamical response of the corona to the heating process 
manifests itself in a comparably slow evolving coronal intensity but rapid changes in Doppler shift. 1 
Subject headings: MHD — stars: coronae — Sun: corona — Sun: UV radiation 


1. INTRODUCTION 

Since the discovery that the outer atmosphere of the 
Sun, the corona, is much h otter than the photosphere 
flCrotriaiill 19391 lEdlenl Il94^) it remains a mystery what 
heats the coronae of the Sun and other cool stars to more 
than a million K. During the last decade a wealth of coro¬ 
nal data was collected through remote sensing, especially 
with the help of the Yohkoh, SOHO and TRACE space 
missions. Since the beginning of coronal physics there 
have been countless suggestions of coronal heating pro¬ 
cesses, as to be found i n the p roceedin gs of conferences 
on coronal heating, e.g.lUlmschneider et al.l 1 1991 1 or the 
SOHO 15 meeting llWalsh et al.ll2004l) . There was never 
a lack of suggestions — the problem is how to prove or 
disprove a model with the help of observations. 

Observations of the corona provide us with the flux 
and energy of photons. To test a model one might use 
an inversion of the observations to derive e.g. temper¬ 
atures and densities, and compare this to the modeled 
plasma parameters. In such an inversion procedure one 
has to make implicit assumpti ons and often the inver- 
sion problem is ill-posed (e.g. Pudge fc Mclntoshlfl99?l 
fMcIntoshli20r^ . The other major approach is forward 
modeling. Based on a model one derives not only the 
plasma parameters, but also the photon spectra result¬ 
ing from the model atmosphere. Like an inversion, also 
a forward model is based on assumptions, but usually 

1 Movies showing the temporal evolution of the synthe¬ 
sized appearence of the model corona similar to Fig.|3 
and Fig.|2l can be found on the web-site http//www.kis.uni- 
freiburg.de/~peter/movie / corona_spec /. 


the assumptions are more easy to control in a forward 
model approach. In the past numerous forward mod¬ 
els have been applied to study various cor onal ph enom- 
ena, like transition region redshifts (e.g. [Hansteen|l993jL 
the c onnectivity of the atmosphere (e.g. |Wikstg^^l 
1998), or catastrophic cooling in loops (e.g. iMuller^^R 
2004 1. Likewise the inversion approach can allow a de- 
tailed comparison, when treated with care. For example 
[Priest et al. i!998l) compared the temperature profiles 
from various loop heating models and compared these 
to the temperature profile derived for a large isolated 
loop as inverted from Yohkoh soft X-ray observations. 

There are many observational challenges a coronal 
model has to explain, and we will here only name a few 
which are of relevance for our work. The structures are 
much more smooth in the corona than in the transition 
region (e.g. iR.eeve Mm. The emission measure, i.e. 
the intrinsic capability for a volume at a given temper¬ 
ature to radiate, is strongly increasing fro m the transi¬ 
tion r egion down to the chromosphere fe.g. lDowdv et all 
mu. The middle transitio n region shows persistent red¬ 
shifts (iDoschek et i _alj ll97fil). while the low corona shows 
a net blueshift~ ifPeteiffi)9fI) . The unresolved motions are 
largest in the middle transition region (e.g. iChae et alJ 
1998af). The temporal variability of the lin e intensities is 
largest in the middle transition region fe.g. lBrkovic et all 
(2003). Here we cannot give a full list of references rel¬ 
evant for coronal heating based on data from Yohkoh, 
SOHO and TRACE, but we would like to emphasize 
the usefulness of spectroscopic investiga tions us in g the 
SUMER EUV spectrometer on SOHO (IWilhelm et alJ 
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IS® 

In the light of the new observations provided in the re¬ 
cent years it is especially encouraging that complex three- 
dimensional ab-initio models for the corona based on 
magneto hydrody namics (MH D) are becoming available 
dGudiksen fc Nordhmdll200^ l2005alb( l. In this work we 
calculate the EUV spectra from these models, which al¬ 
lows us to perform a detailed comparison to observations. 
The aim of this paper is to describe the validity and 
applicability of the approach chosen here to synthesize 
the EUV spectra and to show the huge potential of this 
forward modeling approach. First results o n the spec¬ 
tral sy nthesis have already been published in lPeter et all 
121)1) IK 

We will first give a very quick introduction to the coro¬ 
nal MHD model (Sect.0, before we describe the synthe¬ 
sis of the EUV spectra in Sect.0 The validity and appli¬ 
cability of our approach will be discussed in Sect .0] and 
finally in Sect.0 we will investigate the relation of our 
approach to real solar observations, like the morphology, 
average line shifts and widths or the differential emission 
measure. 

2. 3D MHD MODEL OF CORONAL STRUCTURES 

Here we will only give a very brief overview on the 
basics of the 3D MHD model underlying our spectral 
synthesis. A detailed descr i ption of t he model can be 
found in lGudiksen fc Nordlundl (l2005albH . 

In this model the heating of the corona is due to 
braiding of magnetic flux as a con sequence of p hoto- 
spheric motions, as suggested first by [Parker ( 197 21. The 
braided magnetic field gives rise to currents which are 
then dissipated through Ohmic heating. The horizontal 
motions in the photosph ere are constructed based on a 
Voronoi-tessellation (e.s. lOkabe et allll992 ‘). This flow 
field reproduces the geometrical pattern and the ampli¬ 
tude power spectrum of the velocity and the vorticity. 

The computational box contains a volume of 60 x 
60 Mm 2 in the horizontal directions and 37 Mm vertically 
covering the whole atmosphere from the photosphere to 
the corona with a non-equidistant grid of 150 3 points. 
A 5’th order in space, 3rd order in time, fully com¬ 
pressible 3D magneto-hydro-dynamics (MHD) code on 
a staggered non-equidistant mesh is used, includi ng clas- 
sical heat conduction along the magnetic field llSpitzed 
and a cooling function of optically thin radiative 
losses. A Newtonian cooling scheme is used to keep the 
atmosphere near a prescribed temperature profile in the 
photosphere and chromosphere. 

The simulation starts with an initial condition with a 
magnetic field from a potential field extrapolation based 
on an observed magnetogram, scaled down to fit into 
the computational box. The heating through braiding of 
magnetic flux rapidly leads to a corona which is inter¬ 
mittent in space and time, typically reaching tempera¬ 
tures of about 10 6 K. On average the heating is concen¬ 
trated very much at the bottom of the computational 
domain in the chromosphere and drops exponentially 
towards larger heights. The heat flux into the corona 
of about 2000-8000 W/m 2 is compa rable to the energy 
losses derived from observations fe.g. lWithhroe V Noved 
Ena). The amount of heating in this model is only a 
lower limit, as the heating would stay c onstant or in- 
crease with increased spatial resolution dHendrix et all 


1996 ; Galsgaard fc Nordlund il 996i) . and the initial con¬ 
ditions and the driver are constructed in a way as to 
induce a mi nimum of stress in the magn etic field (see 
discussion in lGudiksen fc Nordlundl2005a h 

For the spectroscopic investigations in this paper we 
use only the last ~20 minutes of the whole simulated time 
span of ~50 minutes, to ensure that our results are suffi¬ 
ciently independent of the initial conditions. Thus in this 
paper the time t=0 minutes refers to a time ~30 minutes 
into the MHD simulation. 

3. CALCULATION AND ANALYSIS OF SYNTHETIC 
SPECTRA 

To study the spectral line properties from the MHD 
calculation we first calculated the emissivities of a num¬ 
ber of emission lines from the transition region and low 
corona at each grid point. Using these emissivities and 
the velocities from the MHD calculation we calculate the 
spectrum for each line at each grid point. Finally we 
integrate the spectra along a line-of-sight in order to ob¬ 
tain spatial maps of line intensity, shift and width. Sim¬ 
ilarly we also compute the average spectra from the box 
to study the average emission line properties. Thus the 
spectra synthesized from the MHD calculation are di¬ 
rectly comparable to observations of EUV spectra, on a 
statistical basis, of course. 

The lines used in this paper are listed in Tabled All 
these lines are optically thin, at least at disk center 
only this allows the concept of synthesizing the spectra 
as outlined above without accounting for radiative trans¬ 
port (see also Sect. 0. The fines have been chosen using 
two criteria. Firstly we wanted to span the whole range 
of temperatures in the transition region from the chro¬ 
mosphere to the corona, i.e. from 10 4 to 10 6 K. And 
secondly the lines should be observable with SUMER to 
allow for a comparison with the observed properties of 
the line profile. 

3.1. Line emissivity 

The line emissivities at each grid point are calcu¬ 
lated under the assumption of ionization equilibrium (cf. 
Sect. l4.2l) . To calculate the various populations of ex¬ 
citation and ionization states we use th e atomic data 
package CHIANT I (Version 4.02; see iDere et al.lll99^: 
lYoung etal] 12001 . Because of computational time we 
first create a look-up table for a large number of tem¬ 
peratures and densities, and then read from this table to 
extract the line emissivities. 

The lines considered in this paper (Tablc0 are pre¬ 
dominantly excited by electron collisions and de-excited 
by spontaneous emission. The emissivity (integrated 
over the line profile, i.e. energy per time and volume) 
is then given by 

e = hvn 2 A 2 i , (1) 

where hv is the energy of the transition, n 2 is the popu¬ 
lation number density of the upper level and A 2 i is the 
Einstein coefficient for spontaneous emission to the lower 
level. Eq. 0 can be rewritten as 

e = G{T,n e ) n 2 e , (2) 

with the electron number density n e and a function 
G'(T, n e ) defined as 

!~urr \ l A n2 n ion Uel H-H , 

G(T , n e ) = hv A 2 1 --— . (3) 

^ion ^el 
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The terms describe the excitation of the line, the ioniza¬ 
tion of the respective ion, the elemental abundance and 
the degree of ionization of the plasma. The number den¬ 
sity of ions in the respective ionization state is 7ii on , the 
number density of the respective element is n e \ and «h 
denotes the number density of hydrogen. 

The first term of ©, hvA 2 i, is a constant given by 
atomic physics. 

Because the excitation of the lines considered here is 
due to electron collisions, the population of the upper 
level ri 2 /riion is basically proportional to the electron den¬ 
sity. This is true only for density insensitive lines and for 
a certain range of densities. Thus the second term of 0 
depends mainly on temperature, but generally also on 
density. We fully account for the temperature and den¬ 
sity dependence when calculating this term using CHI¬ 
ANTI. 

Under the assumption of ionization equilibrium the 
ionization degree n\ ori /n e \ in © depends only on temper¬ 
ature and we use the values calculated bv lMazzotta et all 
as tabulated in the CHIANTI package. Of course, 
this is a strong assumption in a dynamic atmosphere like 
investigated in the present MHD model. Nevertheless 
in large parts of the atmosphere under investigation here 
this assumption can be considered a good first step. This 
will be discussed separately in Sect . 14.21 

In evaluating Eq. 0 we use constant abundances 
Tiei/rin throughout the computational box and adopt the 
most recent solar photospheric values as tabulated in the 
CHIANTI package. Especially we do not account for the 
change of abundances (by a factor of 2 to 3) that is to 
be expected in the chro mosphere, i.e. the FIP effect (e.g. 
Ivon Steiger et al 1 1200(1) . 

In a fully ionized plasma with hydrogen and 10% He¬ 
lium the value of nn/n e is about 0.8. Here we use CHI¬ 
ANTI to calculate this term depending on the degree of 
ionization, i.e. on temperature. 

From this discussion we see that the contribution func¬ 
tion G(T,n e ) as defined in © depends only weakly on 
density when selecting an appropriate line. Thus it is 
often referred to as G(T). 1 Mainly because of the ioniza¬ 
tion equilibrium the function G(T,n e ) strongly peaks in 
temperature. Based on the ionization equilibrium calcu¬ 
lations a canonical value for the width of the contribution 
function is 0.3 in logT (cf. Sect . 14.111 . For density insen¬ 
sitive lines, i.e. G=G(T), the emissivity © reflects that 
the optically thin radiative losses are proportional to the 
density squared. 

To evaluate the emissivity 0 at each grid point of the 
MHD calculation we finally compute the electron den¬ 
sity n e from the mass density p of the MHD calculation 
using the degree of ionization, again as obtained from 
CHIANTI. By an integration along the line-of-sight we 
can then calculate the energy flux density in a given emis¬ 
sion line out of our computational box, e.g. in units of 
W/m 2 These can be easily converted into line radiances 
as observed e.g. at the location of SOHO. 


3.2. Synthetic spectral profiles 

1 Please note that G(T,n e ) is defined differently by different 
authors; e.g. often the abundance is not included in the definition 
of G. 


TABLE 1 

Emission lines synthesized in this study. 




formation temperature: log T [K] FWHM 



ion.eq. 

excitation 

DEM 

MHD 

A logT 

Sill 

1533 

4.60 

4.36 

4.08 

4.25 

- 

Si IV 

1394 

4.81 

4.85 

4.82 

4.90 

0.29 

CII 

1335 

4.67 

4.57 

4.23 

4.64 

0.22 

cm 

977 

4.78 

- 

4.71 

4.84 

0.29 

CIV 

1548 

5.00 

5.00 

5.02 

5.11 

0.25 

OIV 

1401 

5.27 

5.14 

5.15 

5.18 

0.32 

O V 

630 

5.38 

5.35 

5.40 

5.44 

0.28 

O VI 

1032 

5.45 

5.44 

5.50 

5.60 

0.23 

Ne VIII 

770 

5.81 

5.76 

5.82 

5.89 

0.16 

Mg X 

625 

6.04 

6.01 

6.01 

6.06 

0.17 


The line formation temperatures given here are as following from 
different methods: maximum of ionization fraction (“ion.eq.”), 
maximum of ionization fraction accounting for collisional excitation 
(“excitation”), constant pressure DEM inversion using CHIANTI 
(“DEM”) and based on emissivities of the present MHD model 
(“MHD”). The rightmost column gives the width of the respective 
contribution function in temperature as f ollowing from the present 
work in a logarithmic scale. See Sect. I4.1l for a detailed discussion. 


For convenience we calculate all spectral profiles with 
the wavelength given in units of (Doppler) velocities. 

To calculate the emission line profiles at each grid point 
we assume that the line profile has a thermal width with 
according to the temperature T as obtained from the 
MHD calculation at the respective grid point, 


Wth 


2 fc B T \ 1/2 

mi ) 


( 4 ) 


Here fce is Boltzmann’s constant and rrii is the atomic 
mass of the respective ion. 

Furthermore the line is Doppler shifted by the line-of- 
sight velocity at the respective grid point. In this paper 
we concentrate on a vertical line-of-sight, as we would like 
to compare our results of average synthesized Doppler 
shifts to observations of solar disk center Doppler shifts. 
Thus we use the vertical component of the velocity v v 
from the MHD calculation. 

Now the line profile I v at each grid point is given by 


Iv — Ip ea k exp 


(v ~ Vv) 2 \ 

W th ) 


( 5 ) 


For such a Gaussian with an exponential width wth the 
peak intensity / pea k is related to the total intensity Jtot, 
i.e. integrated over wavelength (or here velocity), by 


-ftot — dpeak Wth • (6) 

The total intensity / to t is given through the emissivity 
© as discussed in Sect. rm 
Thus we can use © to calculate the spectrum at each 
grid point as viewed along the (vertical) line-of-sight. To 
obtain the total spectrum we integrate along the line-of- 
sight, 

I° ynth = [ Iv d l . (7) 

J line—of—sight 

This yields two dimensional maps of spectra as they 
would be obtained if an EUV spectrometer like SUMER 
would perform a raster scan. Thus we name these spectra 
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/® ynth , synthetic spectra. These are (in general) not Gaus¬ 
sian, even though in the present study they are mostly 
close to Gaussian. 

The analysis of the synthetic spectra is similar to the 
procedure for observed spectra: we calculate line inten¬ 
sity, position and width of /® ynth at each location of the 
spatial maps. For observed solar spectra it is preferable 
(and almost always necessary) to perform a Gaussian 
fit to derive these parameters, which is mainly because 
of the noise in the data. As our synthetic spectra are 
noise-free and very close to Gaussians it is sufficient to 
calculate the profile moments with respect to wavelength 
(or here velocity v). Then line intensity, shift and width 
are given by 

jsynth = J jsynth ^ ^ (g) 

= ( 9 ) 

*»' yn,h = T^h ( 2 /(•’ - »' yn ‘ t ) 2 C n,h d») (10) 

It is easy to check that this indeed returns the desired pa¬ 
rameters if 7® ynth whould be exactly Gaussian. Examples 
of maps of synthesized Doppler shifts and intensities are 
shown in Fig. 0 (panels a /a and b //3, respectively; panels 
c /7 and d/5 show the box as seen from the sides in inten¬ 
sity; Latin and Greek letters label maps about 7 minutes 
apart in time). 

When analyzing the line width in observations one usu¬ 
ally subtracts the thermal line with as defined in 0 as¬ 
suming the line is formed at a single temperature Tii ne . 
Here we use the temperatures as derived from the analy¬ 
sis of line formation temperatures (see Sect. 14.ll column 
“MHD” in Tabled . This results in the non-thermal line 
width or non-thermal velocity w*( nth characterizing the 
unresolved motions. This includes the velocity fluctua¬ 
tion along the line-of-sight. 

1/2 

< ynth = ^( W synth ) 2 - 2fcB w T ' ine ) • (11) 

It is important to note that when synthesizing the spec¬ 
trum at a given grid point we use the temperature at 
that grid point to evaluate the thermal width following 
0 - After integration along the line-of-sight we subtract 
the thermal width in CD as to be expected at an aver¬ 
age Tii ne in the line formation region. Thus our analysis 
of the data follows closely the analysis of “real” observa¬ 
tional data. 

3.3. Differential emission measure analysis 

The concept of the (differential) emission measure 
(DEM or EM) is widely used in the analysis of stellar 
and solar optically thin plasmas. Thus we will apply this 
tool also to our synthetic spectra. 

Following the discussion in Sect. rrn the line flux (i.e. 
the energy flux of the photons, e.g. in W/m 2 ) at a given 
location on the Sun is given by a height integration of 
the emissivity ©. 

F = J G(T,n e ) n 2 e dh . (12) 


To perform the DEM analysis one has to make the as¬ 
sumption that the integration over height can be substi¬ 
tuted by an integration over temperature, i.e. 

F = J G(T,n e ) DEM dT . (13) 

with the differential emission measure 

DEM = nl ^ . (14) 

It is important to note that this implicitly assumes that 
the temperature varies monotonically with height. The 
highly structured nature of the coro na, o f course, prov es 
this assumption untenable (see Sect. 1X31 and Sect. E3. 

Using a suitable set of lines covering a large range of 
temperatures one can derive the DEM from the observed 
or synthesized line intensities. To account for the density 
dependence of G one has to make an assumption on the 
density in the atmosphere. We use the approximation of 
constant pressure, which is most common for the transi¬ 
tion region because of its small thickness compared to the 
pressure scale height. To perform the DEM inversion we 
employed the CHIANTI package. This also returns the 
respective line formation temperatures (see Sect. l4.ll> . 

We would like to stress here th at the inversion of 
DEM(T) is an ill-pose d problem (e.g. I.Iudge fc Mclntoshl 
Il999t I McIritos hJ 1 200 ()T> and that one has to be very care¬ 
ful in the selectio n of lines, e.g. with res pect to the iso¬ 
electric sequence (iDel Zanna et al.ll20 02). However, our 
emphasis in this paper is not on a best possible DEM 
inversion — we are applying the DEM inversion to our 
data only to show that our synthetic spectra yield a DEM 
distribution with temperature similar to observed ones 
when applying standard techniques. 

4. VALIDATION OF THE APPROACH 

The spectral synthesis as outlined above presents us 
with two major limitations. The formation of lines at low 
temperatures is not well representsed and the ionization 
equilibrium is a strong assumption. 

The MHD model is not including a self-consistent 
treatment of the chromosphere including radiative trans¬ 
port, but is using a Newton cooling mechanism to keep 
the lower part of the computational box at a prescribed 
chromospheric temperature profile. Thus any diagnos¬ 
tics of the plasma at temperatures below, say, 10 4 K is 
meaningless, and will not be touched upon in this paper. 

For the lines formed in the low transition region, below 
logT«4.5, we will be confronted with another problem. 
Because of the very steep rise of density from the tran¬ 
sition region down into the chromosphere, through m 
the formation of these lines is shifted towards consider¬ 
ably lower temperatures, as will be discussed in Sect . 14.II 
Whether this is an artefact of the MHD model used here, 
or a general property of outer stellar atmospheres can 
only be clarified when more complex coronal models be¬ 
come available. If indeed these low temperature lines 
are formed well below the temperature of peak ion frac¬ 
tion (cf. Fig.0i), this would have serious impact on the 
interpretation of observations of these lines. 

For the present study it seems that ionization equi¬ 
librium is a good assumption as outlined in Sect . 14.21 
Despite velocities of some 10 or 20km/s in the mid¬ 
dle transition region this assumption holds. Previous 












Forward modeling of stellar coronae 


5 


studies found contrary results, but they neglected the 
back-reaction of the flow on the temperature gradi¬ 
ent. Future studies with more violent flows, however, 
will have to consider departures from ionization equi¬ 
librium, as is known e.g. from th e modeling of explo¬ 
sive events or reconnecti on (e.g. tRoussev et S3 ISM; 
iRoussev & Galsgaardl il2003) . Those models, however, 
where only 2D-MHD models and the non-equilibrium 
ionization was solved properly only in ID along a line- 
of-sight aligned with the outflow from a reconnection 
site. ID codes with non- equil ibrium ionization exist since 
quite a while (e.g. lHansteen| [1993~). but to move on to a 
full 3D treatment is a major step. 

4.1. Line formation temperatures 

The function G as used in 0 is a function sharply 
peaked in temperature. Thus one can define a line forma¬ 
tion temperature. As G(T ) is mainly determined by the 
fractional density of the ionization stage, one often uses 
the peak temperature of the ionization fraction ni oa /n e \. 
In Table0 values as f ollowing from t he ionizatio n equ i- 
librium calculation of iLandini Kr. Monsignori Foss i. {1199(f ) 
are listed (column “ion.eq.”). One can go one step fur¬ 
ther and also include the temperature dependence of the 
collisional excitation process by using the peak of the 
upper level population in a simple density independent 
wav n ion Mpi • T -1 / 2 exp (—hv/k T) — value s obtained by 
iChae et aIT(ll998al) and lChae et alJ lH998bf) are listed in 
the column “excitation” of Tablc0 However, G depends 
also on density, which is not accounted for by these two 
approaches. A DEM inversion as discussed in Sect . 13.31 
returns the line formation temperatures for a constant 
pressure atmosphere, implicitly assuming a simple ID 
stratified atmosphere (column “DEM” in Tablc0. 

For the analysis of the synthetic spectra from our MHD 
model we do not have to make all these simplifying as¬ 
sumptions, but we can calculate directly the contribution 
function of the respective lines. Employing the proce¬ 
dure outlined in Sect. rm the emissivity is calculated at 
each grid point. Based hereupon we integrate the to¬ 
tal emission from the computational box in each line in 
small temperature intervals (actually, we use intervals in 
log 10 of temperature). As the atmosphere in our model is 
highly structured, the volume each temperature interval 
covers is not necessarily spatially connected. From this 
we can compute the contribution of the respective line 
as a function of temperature. 

These contribution functions are plotted as crosses in 
Fig-Ill for a number of lines. We have used intervals of 
0.02 in log 10 T, but the results change only slightly for 
other values. Except for Si II the contribution in logT is 
roughly represented by a Gaussian. To account also for 
the asymmetries in the contribution function, which are 
due to the highly structured nature of the atmosphere, 
we use the first moment of the contribution function in 
log T to describe the line formation temperature (this is 
practically identical to the center of the respective Gaus¬ 
sian fit, cf. Fig.0. 

A canonical value for the range of temperatures con¬ 
tributing to transition region lines is about 0.3 in logT, 
i.e. 0.3 dex. The full width at half maximum (FWHM) 
of the Gaussian fits is usually a bit smaller than that, 
but 0.3 dex seems to be a good choice, even for a highly 
structured atmosphere as considered here (cf. Fig0and 



log temperature [K] 

Fig. 1.— Contribution to the emissivity as a function of tem¬ 
perature for a number of lines (crosses). Except for Si II the con¬ 
tribution in logT is roughly represented by a Gaussian. Gaussian 
fits are over-plotted as solid lines and the respective full widths 
at half maximum (FWHM) of the Gaussians are given as bars. 
The FWHM values are also given on a log 10 scale with the bars. 
The vertical lines represent the first moment of the contribution 
func tion, used as a line formation temperature in this paper. See 
Sect. l4.l1 


rightmost column of Tablc0 . 

The contribution function for the “coolest ion” dis¬ 
cussed in this paper, Si II, is far from being Gaussian 
or even having a well pronounced peak in temperature 
(Fig.0). There is a small secondary peak just below 
logT=4.6, where the ion fraction peaks (cf. Table0, but 
below logT=4.4 the contribution of Si II rises and well 
exceeds the value at the temperature of peak ion frac¬ 
tion. This is due to the density dependence of the con¬ 
tribution function G (also the DEM inversion using CHI¬ 
ANTI gives a low temperature of log T=4.08 for this ion, 
cf. Table0. The bulk part of the Si II emission origi¬ 
nates from temperatures below logT=4.4, where the ion 
fraction of Si + is still low. As there is no self-consistent 
treatment of the upper chromosphere, i.e. for temper¬ 
atures up to some 20 000 K (logT=4.3), the diagnostic 
potential of Si II is limited within the framework of this 
study. 

The main point of this discussion on Si II is to show 
that based on the present model, it seems not very use¬ 
ful at all to use Si II or other comparably cool ions to 
extract information on the state of the transition region. 
The temperature regime below 20 000 K, where also Ly- 
a is formed, is not well described by a simple stratified 
optically thin atmosphere as assumed for (almost all) ob¬ 
servational investigations of EUV emission lines. 

4.2. Ionization balance 

To properly account for ionization effects one would 
have to solve the rate equations for each species, i.e. for 
a considerable number of ionization and excitation states. 
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In general for an ionization/excitation state i this is 

+ V • {m v) = {-la + 7 ji) , (15) 

3 

where rii is the number density and v the bulk velocity. 
The summation on the right-hand-side is over all loss 
processes to states j and gain processes from j with rates 
7 ij and 7 ^. 

Currently we cannot solve these rate equations parallel 
to the 3D MHD problem, simply because of computa¬ 
tional time. 

Therefore, when calculating the emissivity of the var¬ 
ious lines as outlined in Sect. rm we assume ionization 
equilibrium. To check this severe simplification we com¬ 
pare the ionization times to dynamic times given by the 
temperature gradient and the velocity of the plasma. 

4.2.1. Ionization times 

To calculate the ionization times 7 i on at each grid point 
in the box of the MHD calculation we use the ionization 
rates Cion (P) as parameterized in llArnaud & Rothenflud 
CMS)- These rates are still up-to-date and, except for 
Fe, are also used in the recent ioniz ation equilibrium cal¬ 
culations of iMazzotta et all 1. 1998 1). We account for di¬ 
rect ionization as well as for excitation-autoionization, 
where appropriate. 

At each grid point of the MHD model we first evaluate 
the ionization rate Ci on for the respective temperature 
and obtain the ionization time by multiplying with the 
electron density from the MHD model, 

- = Tie Cion (16) 

Hon 

This is done for each ion of the lines listed in Tabled 

As we are interested in the effects of the (non¬ 
equilibrium) ionization on the emission lines, we calcu¬ 
late an ionization time for each line by considering only 
those grid points of the MHD calculation with tempera¬ 
tures ± 0.1 dex in logT from the line formation tempera¬ 
ture Tii ne . We use values of Ti; ne as derived in the previ¬ 
ous section (cf. column “MHD” in Tabled . The range of 
temperatures was chosen to be ±0.1 dex, i.e. 0.2 in logT 
according to the widths of the contribution function (cf. 
FWHM in Tabled. 

This leads to a distribution of ionization times in the 
line forming region of the respective line. For the case of 
C IV (1548 A) this histogram is shown as a dashed line 
in Fig. EH As the ionization time for this line we use 
the median value of the distribution (solid diamond). To 
characterize the width of the distribution we compute 
the standard deviation of the ionization times on a log¬ 
arithmic scale. This is shown as a bar in Fig. EH. The 
ionization times for the other lines are displayed as just 
those bars with diamonds for the median values in Fig.El 

4.2.2. Dynamic times 

We use two procedures to determine a dynamic time, 
the first closer linked to the picture of particles getting 
ionized while flowing up a temperature gradient, the sec¬ 
ond more formally bound to the advective term in the 
rate equations. 

The flow time Tfl ow we define as the time needed for a 
test particle with velocity v to cross a given temperature 



Fig. 2.— The left panel (a) shows a histogram of ionization 
times Ti on for C 3+ (dashed) as compared to histograms of flow 
times Tfl ow (solid) and advective times T a dv (dashed-dotted) at the 
temperatures where C IV (1548 A) is formed. The diamonds show 
the median values of the respective times and the bars indicate the 
standard deviations. 

The right panel (b) shows the histograms of the ratios of flow to 
ionization time, T'low ,/ r ion (solid), an d ad vective to ionization time, 
Tadv/tlon (dashed-dotted). See Sect. IOI 



log 10 temperature [K] 


Fig. 3. — Comparison of time scales for ionization, flow and 
advection. The diamonds show the median values of the ionization 
time Tion of the respective lines (standard deviation shown as bars, 
cf. Fig.|§Ji). The solid line shows the median flow time Tfl ow at 
the respective temperature in the MHD model. The dotted lines 
indicate the standard deviation of Tfl ow . In addition the dashed- 
dotted li ne s hows the median values of the advective times r a d v - 
See Sect. 101 


difference (AlogT)i on . Thus within the time Tfl ow we 
put the test particle to a higher temperature and ask if 
the particle is ionized after this time. If this is not the 
case, i.e. if Tfl ow > Ti on , ionization equilibrium is a good 
approximation. 

We choose a value of (A log T) j on = 0.1, as this is about 
the half width at half maximum of the contribution func¬ 
tion for the emission lines (FWHM/2 in TablcEJ. When 
moving an ion by such a temperature difference it just 
starts “seeing” another ionization equilibrium. 

Thus the inverse of the flow time is given by 


1 V(logT)--f) 

Tfl ow (A logT) 

ion 


(17) 


where the temperature gradient and the velocity vector 
v are taken from the MHD model, v is the unit vector 
along v. The temperature gradient is projected on the 
flow direction as we are interested only in changes along 
the path of the test particle. Furthermore we consider 
only locations where the flow is in the direction of in- 
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creasing temperature as only this will be affected by the 
ionization processes. 

Like for the ionization times we now compute the dis¬ 
tribution of flow times in each line formation region. 
These histograms can be directly compared to those of 
the ionization times, as it is done in Fig.|2li for C IV. 
This shows that for C IV the flow time is larger than the 
ionization time, which is also confirmed by the histogram 
of ratios of flow to ionization times (Fig.®)). 

In Fig.®the median values of the flow times, Tfl ow , are 
shown also for other temperatures in the MHD model 
as a solid line. The dotted lines indicate the standard 
deviation of Tfl ow . As expected the flow time is small¬ 
est in the middle transition region below 10 5 K as there 
the temperature gradient is steepest. Tfl ow increases by 
almost two orders of magnitude towards the corona as 
well as to the chromosphere. One should note that the 
absolute values of the flow times depend on the choice of 
(A log T) ion . 

A more direct way to compute a dynamic time is to 
evaluate the left-hand-side of the rate equations 11151) . 

If considering only the process of electron collisional 
ionization for a single ionization state, the rate equation 
ltTT)l) reads 

- + - V • (nr) = —n e C ion , (18) 

n df n 


with rii = n being the density of the respective ion, rij = 
n e the electron density and 7 = rii n e Cion the ionization 
rate. The right-hand side is the inverse ionization time 
as used in da and from the left-hand-side we can define 
an advective time scale 


— = - |V • (nu)| , 

Tadv 


(19) 


The divergence of the particle flux density V-(nv) is cal¬ 
culated from the MHD model. Just like for Tfl ow we have 
computed the median values of r a d v and over-plotted 
them in Fig.® as a dashed-dotted line. At all temper¬ 
atures the advective time scale is at least as large as the 
flow time scale defined above, r a d v > Tflow- 


4.2.3. Comparing ionization and dynamic times 

A comparison of the ionization time r; on with the dy¬ 
namic times Tflow and r a( jv shows that in general the 
assumption of ionization equilibrium is satisfied in the 
model under investigation here. 

A more detailed comparison of the distribution of ion¬ 
ization and dynamic times shows that there are regions 
where one definitely should not use ionization equilib¬ 
rium. Mostly, like for C IV only a small fraction of the 
volume has flow times smaller than ionisation times (cf. 
in Fig.®> only a minority has values of Tfl ow /rio n <l)- 
In our study it is only for C III and Si IV that the as¬ 
sumption of ionisation equilibrium is violated in the sense 
that Tfl ow <Tion> but even then the ionisation time is still 
smaller (or comparable) to the the advection time T a( j v 
(cf. Fig.®. 

From this we can conclude that for the present work 
and as a first step for a spectroscopic analysis of complex 
3D MHD coronal models we might well use ionization 
equilibrium. Future work, however, should try to include 
also non-equilibrium effects. 



Fig. 4.— Average temperature gradient along the streamlines 
in the MHD model underlying the present analysis as a function 
of temperature for regions with flow speeds of about 10 km/s 
(dashed) and 0.5 k m/s (so li d). F or comparison the gradient in 
the static model of I Gabrieli 01976D is shown as a dot-dashed line. 
The bar repr esents the scatter of gradients in each temperature 
bin. See Sect. 14.2.41 


4.2.4. Ionization equilibrium and flows 

It is often argued that departures from the ionization 
equilibrium become of vital importance as soon as the 
velocities across a temperature gradient become “large 
enough”. As shown in the preceding subsection in our 
computation the ionization times are (generally) still 
shorter than the time scales of the flows, and hence the 
assumption of ionization equilibrium is not too bad. 

T his is in c ontras t to earlier considerations, e.g. of 
Uoselvn et, all ( 1979 1. who argued that velocities of the 
order of 10 km/s will certainly lead to a violation of ion¬ 
ization equilibrium. They used the t emperature gradi¬ 
ents from static models (e.g. lGabrielll976l) to calculate 
time scales of the plasma flowing across these tempera¬ 
ture gradients. However, such an analysis does not ac¬ 
count for the back reaction of the flow on the temperature 
gradient, and thus a too small time scale is derived. 

To illustrate this we plot the average temperature gra¬ 
dient (along streamlines) as a function of temperature 
within the box of the MHD model, once for only those 
regions with velocities of about 10 km/s (dashed) and 
about 0.5 km/s (solid) in Fig.® It is clearly evident that 
the temperature gradient is smaller in regions with higher 
velocities. For comparison we plot the temperature gra¬ 
dient in a static co ronal (funnel) model (dot-dashed in 
Fie.m lGabrielU97filh In the case of small velocities the 
temperature gradient in the MHD model underlying the 
present analysis is almost as steep as in a static model 
(factor ~2 in the middle transition region). As there is 
no heating in the transition region of the static model 
of iGabriell 119701) . while in the 3D MHD model with 
flux-braiding the heating is strongly concentrated at low 
heights, one expects this somewhat shallower gradient in 
the latter case. The main point to stress here is that if the 
flow speed is larger, the gradient becomes much flatter, 
by about an order of magnitude, as the flow is smooth¬ 
ing the temperature gradient (cf. solid and dashed line 
in Fig.®). 

From this discussion one might draw the following con¬ 
clusion. For a higher velocity the temperature gradient 
will decrease, which partly compensates the effect of a 
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faster speed when calculating the flow time using e.g. 
EJ. By this the back reaction of the flow on the tem¬ 
perature helps in simplifying the treatment of ionization 
to the case of equilibrium. However, for models with 
more violent flows, the assumption of ionization equilib¬ 
rium certainly will break down. 


4.3. DEM from the MHD calculations 


To check the validity of the inversion of the differen¬ 
tial emission measure (DEM) discussed in Sect . 13.31 we 
calculate the DEM directly from the MHD results, i.e. 
the density. For a traditional DEM inversion one has 
to assume a ID atmosphere with a monotonic increase 
in temperature. Only then the DEM definition in o 
makes sense. However, this is not the case in a highly 
structured atmosphere as under investigation here. Thus 
we start by defining a volume emission measure 

EM V (T) = f n 2 e dV , (20) 

Jv 

where the integration is over a volume typically con¬ 
tributing to an emission line formed at temperature 
T. According to the contribution functions to the line 
CSect. nTl) the volume covers a temperature interval of 
±0.15 in logT, i.e. AlogT=0.3. 

To move to a height-related expression the volumet¬ 
ric emission measure is divided by the area A of the 
box giving an average column emission measure EM/, = 
EMy/A. Now one can formally substitute the volume 
integration by a height integration d h = dV/A and this 
by an integration over temperature, 

EM/, = J nl dh EM/, = J n\ ^ dT . ( 21 ) 

Using the definition of the DEM o one gets 


EM/, 


DEM dT . 


( 22 ) 


Thus one can approximate the DEM by 


DEM M hd 


EM,, 

AT 


EMy 
A AT ’ 


(23) 


where AT is the temperature range corresponding to 
AlogT=0.3, i.e. AT increases with temperature. 

Using PH this allows the calculation of the DEM di¬ 
rectly from the MHD model, which can be compared to 
the DEM as resulting from the inversion of line intensi¬ 
ties as outlined in Sect. lT3l 

In Fig. 0 we plot the DEM derived directly from the 
MHD model as following PH as a dot-dashed line along 
with the DEM from the inversion using the line emissiv- 
ities (solid) as discussed in Sect. 13.31 We see that they 
compare relatively well, even though the difference at 
high temperatures is noticeable. As this discrepancy is 
probably (at least partly) due to the highly structured 
nature of the atmosphere, a further discussion is shifted 
to Sect . 15.5| 


5. RESULTS 

Here we will only give some examples of the results 
that may be achieved using the spectroscopic analysis 
described in this paper. More specific problems, such as 
time variability, evolution or spatial structure will be ad¬ 
dressed in future investigations and will provide new tests 
for the model of coronal heating through flux braiding. 


5.1. Morphology: intensity and Doppler shift maps 

By investigating maps in intensity and Doppler shift 
in various lines when integrating along a line-of-sight 
through the box, one can now study the morphology of 
the transition region and corona. Such maps are pro¬ 
vided in FigP for C IV formed in the low transition 
region at ~10 5 K (top rows labeled with 1) and for Mg X 
formed in the low corona at ~ 10 6 K (bottom rows la¬ 
beled with 2). There the respective leftmost panels (a/a) 
show the Doppler maps when looking from straight above 
and the panels b//3 show the same in intensity — this 
would correspond to an observation at disk center. The 
right panels (c/ 7 , d/<5) show intensity maps when looking 
from the sides at the box (displayed with correct hori¬ 
zontal/vertical aspect ratio) — this corresponds to the 
appearance at the limb. For each line the snapshots at 
two times are shown at some 4 and 11 minutes into the 
simulation (labeled by Latin and Greek letters, respec¬ 
tively). 

In an earlier publication we reported shortly on the 
more structured appearance of the transition region be¬ 
ing due to heat conductio n an d t he n u mero us structures 
in coronal Doppler shifts llPeter et al.ll2004 1. 

When comparing the intensity maps of the two time 
steps only seven minutes apart, it is striking how much 
the transition region has changed in this time, while the 
corona still looks pretty much the same in intensity. This 
is of special interest, if one recalls that the driver of the 
corona, the photospheric granular motion, acts on the 
time scale of some 5 to 10 minutes. Thus the transi¬ 
tion region directly responds to the changes in the pho¬ 
tospheric magnetic structure, while the corona is reacting 
much slower. This difference is mainly due to the large 
difference in cooling time scale and the very efficient re¬ 
distribution of energy trough heat conduction at high 
temperatures. 

Not surprisingly also the Doppler map of the transi¬ 
tion region changes significantly during the only seven 
minutes on small scales. Despite the smooth, slowly- 
evolving appearance in intensity, the corona reveals its 
dynamic structure in the Doppler map. While the heat 
conduction quickly distributes the deposited energy, still 
the dynamic response of the plasma to the heating pro¬ 
cess is present, showing up as flows causing the Doppler 
shifts. And just like in the transition region, the coro¬ 
nal Doppler shifts show large variations on small spatial 
scales within the time scale of the photospheric driver. 

This high variability of the coronal Doppler maps 
shows the need to study coronal variations on short time 
scales not only using intensity maps like provided by 
EIT/SOHO or TRACE. If one could obtain Doppler 
maps with high temporal resolution, one would have 
access to the dynamic response of the corona to the 
heating process, which is not available through intensity 
maps alone. Unfortunately current slit spectrographs 
like SUMER or CDS on SOHO do not provide maps of 
sufficient size in a reasonable time cadence. 

5.2. Vertical structure: roiling and loops 

To show the further potential of the forward modeling 
technique presented here, we investigate a vertical slice 
of the computational box cutting the two main polarities 
forming the active region. In Fig. ED c, d, and e we show 
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Fig. 5.— Spatial maps in Doppler shift and intensity in the lines of C IV (1548 A) and Mg X (625 A) formed at about 10 5 and 10 6 K, 
respectively, for two different time steps of the MHD model about 7 minutes apart (at t =4 min and 11 min). 

The two top panels labeled with 1 show C IV, the two bottom panels labeled with 2 show Mg X. The maps of the earlier time step are 
labeled with Latin letters (1st and 3rd row), those for the time step 7 minutes later with Greek letters (2nd and 4th row). 

The panels (a/a) show the Doppler shift of the synthesized spectra as seen from straight above, the panels (b//3) show the same for the 
intensity of the respective line. This corresponds to the appearance near disk center. Panels (c/ 7 ) and (d/<5) show side views of the 
computational box along the x and y axis in line intensity, which resembles the appearance at the limb. 

The intensities I are scaled with respect to the average (median) intensity (I) of the respective map. 


the location (between the lines) of the vertical slice in 
the horizontal plane for the intensity maps in C IV and 
Mg X as well as the photospheric magnetogram. The left 
panel (Fig.®) shows the vertical slice, the emission from 
C IV and Mg X displayed on top of each other. Please 
note that in contrast to Fig. Inhere the aspect ratio is not 
unity, but about a factor of two. 

The emission from these two lines comes from practi¬ 
cally disjunct regions. The coronal Mg X emission stems 


from the diffuse loop-dominated part in the upper part of 
the image, while the transition region C IV emission orig¬ 
inates from a very thin layer below the corona. Panel b 
of Fig. El shows the area from ~2-8 Mm height vertically 
stretched to show the transition region more clearly. At 
least in this region the transition region emission comes 
from an area below the coronal source region, i.e. from 
the footpoint regions of the coronal magnetic structures. 

Fig® also reveals that the transition region is roil- 
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horizontal coordinate X [Mm] —I i n the photosphere 


Fig. 6 .— Composite view of a thin vertical slice of the computational box in C IV (~10 5 K) and Mg X (~10 6 K) at time £=4min. The 
left panel (a) shows the vertical slice as defined by the black lines in the top views (c,d) when integrating along the Y direction; actually 
panels (c,d) are the same as panels (lb,2b) of Fig.JS] The location of the slice with respect to the photospheric magnetic field map is shown 
in panel (e). The middle panel (b) shows part of panel (a) stretching the vertical axis to highlight the structure in the low parts of the 
atmosphere. The emission as seen in C IV and Mg X comes from disjunct areas of the v ertica l slice. The lower “line” shows the emission 
from C IV, while the upper more diffuse emission stems completely from Mg X. See Sect. o 


ing field, while being very thin in the vertical direction 
(some 100 km). This is to be expected, as the heating 
is highly variable in space and time and thus the loca¬ 
tion of the transition region is to be expected to move up 
and down to adjust the coronal pressure in accordance 
with the he ating rate. This is well known from ID loop 
m odels (e.g. lHansteenj 1 9931 1 and also commented on by 
iGudiksen fc Nordlundl (2005a). Despite being very thin, 
the roiling of the transition region leads to a relatively 
thick appearance when integrating along a horizontal di¬ 
rection. In Fig. EJ> the roiling of the transition region has 
an amplitude in height of about 1 Mm or more, while 
at other places and times it can be even more. Thus 
the roiling leads to a transition region that would appear 
several Mm thick when observed at the limb, as nicely 
illustrated by the side views in Fig.0 panels lc and Id. 
This corresponds well will observations of lim b intensity 
profiles at the limb (e.g. iMariska et alllll978 1. 

5.3. Average line shifts 

The time-average of the line shifts derived from the 
synthe tic corona presen ted here were discussed previ¬ 
ously ijPeter et all 1200411 . In Fig.0 we plot the aver¬ 
age line shifts as a function of line formation temper¬ 
ature for the two individual time steps, seven minutes 
apart, which have been displayed also in the spatial 
maps in Fig.0 The Doppler shifts are calculated for 
a vertical line-of-sight, i.e. when looking at the com¬ 
putational box from straight above. The height of the 
bars and rectangles represents the spatial scatter for each 
time step (1/2 of the standard deviation). The trend 
found in observ ations is plotted as a thick dashed line 
(compiled from iBrekke et all 119971 iChae et all |l998b; 
Peter fe ludgell 199(1. ~ 

It is obvious that the average Doppler shifts change 
quite dramatically, and the overall shape varies from be¬ 
ing very close to the observations (bars) to a more flat 
profile with comparable Doppler shifts at all tempera¬ 
tures (rectangles) — and this happening in only seven 
minutes. The time variation of the Doppler shift for 
some 20 minutes is shown in Fig. llOb for a number of 
lines. Here we see that for all lines, also for the coronal 
Mg X line, the variation over these 20 minutes is compa- 
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Fig. 7.— Comparison of synthesized and observed Doppler shifts 
for two time steps about 7 minutes apart (same times as in Fig.151 
t= 4,11 min). The diamonds show the Doppler shifts of the average 
spectra and the bars indicate the standard deviation of the Doppler 
shifts of a spatial map as seen looking from straight above on the 
computational box for the earlier time step. The rectangles show 
the same but for 7 minutes later. The change in this short time is 
noticeable. The lines are plotted at the formation temperature as 
following from the contribution to the emission from the present 
work (cf. TablcQ. The t hick dashed line shows the trend as found 
in observations. See Sect. c~n 


rable or larger than the spatial scatter in the data. This 
again emphasizes the potential provided by observations 
of Doppler shifts throughout the outer solar atmosphere, 
as pointed out at the end of Sect. IQ 
In the present model the persistent redshifts in the 
tr ansitio n region, puzzling theorists since their discovery 
bv lDoscliek et all ( 1976 1. are caused by the flows induced 
by the heating through braiding of magnetic flux. Being 
free of spurious assumptions, this is the first time one 
finds an explanation of the Doppler shifts providing a 
qualitative and quantitative representation of the red- 
shifts in the middle transition region. The blue shifts 
towards higher temperatures are not matched yet, but 
we see a clear indication that above logT«5.2 the red- 
shifts are decreasing, and at some periods of time in the 
model run, we are quite close to the observations (cf. 
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Fig. 8. — Comparison of synthesized and observed non-thermal 
widt hs. Similar as Fig.Q but now for non-thermal widths. See 

Sect. |^3] 


bars in Fig.0. Further investigations, especially with an 
improved upper boundary condition or a larger box, will 
have to address this discrepancy. 

5.4. Average non-thermal line widths 

Another important test for any coronal model is pro¬ 
vided by a comparison of the non-thermal line widths. 
These contain information on the non-thermal unre¬ 
solved motions, introduced by the limited resolution in 
space and time of any instrument. Averaging in space 
and time will mix spectral information with different ve¬ 
locities and will lead to a broadening of the line. It is 
important to note that also the integration along the 
line-of-sight contributes to the non-thermal broadening, 
and there is (almost) nothing we can do about this from 
an observers point of view. 

In Fig.|Sl we plot the non-thermal line width of the 
average spectrum for the two time steps seven minutes 
apart already displayed in Fig.[S]and[7| (again when look¬ 
ing at the computational box from straight above). The 
thick dash ed line shows the t rend in observati ons com¬ 
piled from lChae et alJ (Il998af) and iPeteil (|2001f) . 

In the corona and the low transition region the 
non-thermal width from the synthesized spectra match 
roughly the observed values, but in the middle transition 
region the synthetic spectra show significantly smaller 
widths than the observations. There is (at least) one pos¬ 
sible explanation for this. The high non-thermal broad¬ 
ening observed on the Sun might be due to small scale 
velocities connected with the heating process itself, i.e. 
the nano-flares induced by the flux-braiding. And this 
effect is to be expected to have the largest effect in the 
middle transition region, where the time scales are short¬ 
est. Because of their tiny spatial scale, these nano-flares 
cannot be modeled in the MHD simulation, which the 
spectral synthesis is based upon. Future models with 
increased spatial resolution will have to show if this mis¬ 
match with the observations can be resolved. 

5.5. What do we learn from the differential emission 
measure? 

As outlined in Sect. Id.dl and rOl we perform an differen¬ 
tial emission measure (DEM) analysis using the spectral 


Fig. 9. — Comparison of the differential emission measure (DEM) 
as derived from the synthesized and observed spectra. The solid 
line shows the fit from the D EM i nversion based on the synthesized 
lines displayed as bars (Sect. ESI for time t=4min). The lines are 
displayed here at the formation temperature as following from the 
DEM inversion (cf. TableQJ. This differs slightly from the tem¬ 
peratures used in Figs. |T]and[H] The rectangles show the situation 
about 7 minutes later at £=11 min, the changes being very small. 
The thick dashed line is based on a DEM inversion using observed 
quiet Sun disk center line radiances observed with SUMER, scaled 
by a factor of two. The dot-dashed line displa ys th e DEM cu rve a s 
following directly from the MHD results (Sect. ehj See Sect. EH 


lines listed in Table^using the atomic data package CHI¬ 
ANTI. The resulting DEM curve for one single time step 
(at 4 minutes into the simulation) is plotted in Fig^as 
a solid line. The diamonds show the DEM at the line 
formation temperature of the respective line, multiplied 
by the ratio of the observed emissivity to the one pre¬ 
dicted by the DEM fit. Thus the displacement of the 
diamonds from the solid line represent the error of the 
DEM fit. The height of the bars represent the scatter of 
the emissivities. 

As we discussed previously (iPeter et all 1200 41 this is 
the first model to reproduce the overall shape of the 
DEM curve quantitatively and qualitatively, especially 
the increase of the DEM towards lower temperatures be¬ 
low logT=5.3. For comparison we plotted an DEM in¬ 
version based on radi ance data at disk center taken from 
iWilhelm et al] dl99Sh . for Mg X and Si II we re-evaluated 
some SUMER disk center spectra. The agreement is re¬ 
markable, except for the highest temperatures. This dis¬ 
crepancy is partly due to a lack of constraining lines at 
higher temperatures, as the maximum temperature in 
the MHD simulation is about logT«6.3. Thus the DEM 
inversion is not very well defined at high temperatures. 
This becomes especially important when discussing the 
temporal variability of the DEM below. 

The DEM inversion based on the synthesized spectra 
roughly agrees with the DEM derived directly from the 
MHD models (cf. Sect. 14.3D . While for the DEM inver¬ 
sion one implicitly assumes a simple ID stratified static 
atmosphere (cf. Sect. l3.3(l this is certainly not the case 
in the computational box. In this light it is even sur¬ 
prising that the DEM derived directly from the MHD 
model (Fig.® dot-dashed line) roughly compares to the 
inversion, especially below logT=5.0. 

In Fig.HJthe rectangles represent the situation concern¬ 
ing the DEM inversion 7 minutes after the time step dis- 
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Fig. 10. — Temporal evolution of the differential emission mea¬ 
sure (DEM) and the average Doppler shift. The top panel (a) 
shows the evolution for the DEM at three different temperatures 
as following from DEM inversion fits like the one shown in Fig.f^l 
The bottom panel (b) shows the variation of the Doppler shifts 
for emission lines spanning temperatures from logT~4.6 to 6.0. 
Please note that the variation in DEM in much smaller than the 
uncertainties (represented by the bars in both panels) and is thus 
not significant. In contrast, the variation in Doppler shift is signif- 
ican t. Th e Figs. -|3 show snapshots at £=4min and 11 min. See 
Sect . 18.51 


cussed so far (represented by the bars). The change is 
very small. This is emphasized by the temporal evolution 
of the DEM inversion fit at temperatures representing 
the low and middle transition region as well as the low 
corona lFig.11 Oh.b It is clearly evident that the changes 
in the DEM fit are only very small. At logT=5.9 the 
fit is “jumping”, as the DEM inversion is not well con¬ 
strained at high temperatures (see above). Thus also the 
jump of the fit for logT=5.9 at time ~18 minutes is an 
artefact and not real. 

The important result here is that the DEM hardly 
changes with time (over the 20 minutes of the simula¬ 
tion), but within the uncertainties stays constant. This 
is in strong contrast to the Doppler shifts, which change 
significantly at all temperatures. 

The DEM curve thus provides an important test for a 
coronal heating model, in particular in the sense that it 
is not easy to get the overall shape right. However, once 
one gets the overall shape right, it seems that the DEM 
stays fixed, and therefore provides only limited potential 
for further diagnostics of the structure and dynamics of 
the coronal heating process. 


6 . CONCLUSIONS 

We have presented a procedure how to synthesize EUV 
spectra from a complex 3D MHD model for a stellar 
corona, in which the heating is due to flux braiding 
through photospheric footpoint motions. An investiga¬ 
tion of the validity of the ionization equilibrium showed 
that in the present case this simplification is justified. 
Because the heating is concentrated very low in the at¬ 
mosphere, the temperature gradients are much shallower 
than in traditional static models, and thus the typical 
dynamic (or flow) times are larger than the ionization 
times. However, for models with more violent flows one 
will have to properly solve the ionization rate equations. 
We found that due to the strong heating at the foot- 
points of the corona and the resulting strong increase of 
density towards the chromosphere, cool lines like from 
Si II tend to be formed well below their ionization equi¬ 
librium temperature. Thus the cool lines are of limited 
diagnostic value for the transition region. 

Together with the spectral synthesis the MHD simu¬ 
lation provides a forward model, which allows a direct 
comparison to observations, and by this is a powerful 
tool to investigate the corona. In this paper we describe 
some of the potentials of this method. 

The spatial maps of intensity show a much higher spa¬ 
tial and temporal variability in the transition region than 
in the corona. In the transition region the temporal vari¬ 
ability occurs on a time scale compatible with the pho¬ 
tospheric driver of the coronal heating. In the Doppler 
shift maps, however, the spatial variability is also very 
large in the corona, much stronger than the variability in 
intensity. This shows the need for further instrumenta¬ 
tion to get access to coronal Doppler maps with sufficient 
temporal resolution (below 5 minutes). 

Likewise the time variation of the Doppler shifts is sig¬ 
nificant, while the differential emission measure (DEM), 
derived from the intensities, shows only a very small vari¬ 
ation. While this model is the first to reproduce qualita¬ 
tively and quantitatively the form of the DEM variation 
with temperature, it also shows the limitations of the 
DEM analysis in a highly dynamic atmosphere. To un¬ 
derstand the dynamics of the corona it is vital to use the 
Doppler shifts as a crucial test for the model, and not 
only the line intensity or emission measures. 

This study shows the pivotal importance of forward 
modeling for our understanding of stellar coronae, as it 
provides numerous controllable ways to compare model 
results and observations. 
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